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CRITICAL PARTITIONS AND NODAL DEFICIENCY OF 
BILLIARD EIGENFUNCTIONS 

GREGORY BERKOLAIKQi, PETER KUCHMENTi, AND UZY SMILANSKY^.s 



PsJ ■ Abstract. The paper addresses the nodal count (i.e., the number of nodal 

domains) for eigenfunctions of Schrodinger operators with Dirichlet bound- 
ary conditions in bounded domains. The classical Sturm theorem states 
that in dimension one, the nodal and eigenfunction counts coincide: the n- 
^^ th eigenfunction partitions the interval into n nodal domains. The Courant 

Nodal Theorem claims that in any dimension, the number of nodal domains 
Vn of the nth eigenfunction cannot exceed n. However, it follows from an 
ly-v . asymptotically stronger upper bound by Pleijel that in dimensions higher 

^SJ ' than 1 the equality can hold for only finitely many eigenfunctions. Thus, 

in most cases a "nodal deficiency" dn = n~ z/„ arises. One can say that the 
nature of the nodal deficiency has not been understood. 
O ■ It was suggested in recent years that, rather than starting with eigen- 

functions, one can look at partitions of the domain into i' sub-domains, 
asking which partitions can correspond to eigenfunctions, and what would 
2 ' be the corresponding deficiency. To this end one defines an "energy" of a 

partition, for example, the maximum of the ground state energies of the 
sub-domains. One notices that if a partition does correspond to an eigen- 
function, then the ground state energies of all the nodal domains are the 
^ ■ same, i.e., it is an equipartition. It was shown in a recent paper by Helf- 

0^ . fer, Hoffmann-Ostenhof and Terracini that (under some natural conditions) 

QO ■ partitions minimizing the energy functional correspond to the "Courant 

sharp" eigenfunctions, i.e. to those with zero nodal deficiency. 

In this paper it is shown that it is beneficial to restrict the domain of 
r_^ , the functional to the equipartition, where it becomes smooth. Then, under 

some genericity conditions, the nodal partitions correspond exactly to the 
critical points of the functional. Moreover, the nodal deficiency turns out 
to be equal to the Morse index at the corresponding critical point. This 
explains, in particular, why the minimal partitions must be Courant sharp. 

pi. 1. Introduction 

We consider the Schrodinger operator 

(1) H = -A + V{x) 

with Dirichlet conditions in a connected bounded domain Q C M.'^. We will 
assume that the domain has a smooth boundary and that the real potential V 
is also smooth on Q. While these assumptions are overly restrictive, we do not 
want to burden our considerations with less significant details. See the final 
section for additional remarks. 

The operator H can be defined via its quadratic form 
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h[u,u]= / \\/u{x)\'^dx + / V{x)\u{x)\'^dx 
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with the domain Hq{Q). Thus defined, it is self-adjoint in L2(fi) and has real 
discrete spectrum of finite multiplicity 

Ai < A2 < A3 < . . . , 

where lim A„ = cxd. It has an orthonormal basis of real- valued eigenfunctions 

ipn such that ipi{x) > 0. We will sometimes use the notations H{Q) and Xj{fi), 
when we need to emphasize the dependence of the operator and its spectrum 
on the domain. 

Definition 1. For a function f{x) we will be interested in its nodal (zero) 
set 



zif) ■■= r\o) = {xen I fix) = 0}. 

The complement Q\N{f) is the union of connected open sub-domains Di, . . . ,D^ 
of Q, which we will call nodal domains. The nodal domains form the nodal 
partition P{f) = {Dj} corresponding to the function. 

We will be mostly interested in the case when f{x) is an eigenfunction ipn, 
and thus in its nodal set Z{ipn), its nodal domains Di, . . . , D^ of f2, and its 
nodal partition P{ipn) = {Dj}- 

A lot of attention has been paid to the nodal structure of eigenfunctions 
(e.g., [T5]4T9t l30 y 43 [ H6] and references therein), and in particular to the number 
z/ (or z/^^, if one wants to emphasize dependence on the eigenfunction) of nodal 
domains of the n-th eigenfunction ipn oi H. 

In spite of more than 300 years history of this topiqj, open questions still 
abound. We will discuss here one of them, the issue of the so called nodal 
deficiency. The classical Sturm theorem states that in dimension one, the nodal 
and eigenfunction counts coincide: z/^^^ = n. The Courant Nodal Theorem [TU 
Vol. I, Sec. V.5, VI. 6] asserts that in any dimension, the upper bound on the 
number of nodal domains still holds: 

While z/^^ = 1 and u^^ = 2, it follows from an asymptotically stronger upper 
bound for u^^ by Pleijel [IH] that in higher dimensions the equality u^^ = n can 
hold only for a finitely many values of n. Moreover, there are known example 
of eigenfunctions ipn with arbitrarily large index n that have just two nodal 
domains. 

The eigenfunctions tpn for which z/^^ = n, are sometimes called Courant 
sharp. The non- negative difference 

dn-=n- z/^„ 

is said to be the nodal deficiency of an eigenfunction. It is believed that the 
integer sequence {dn} contains much information about the geometry of the 
domain (see, for example, [SlIiniEZlEZlES])- However, we are not aware of any 



Robert Hooke observed on 8 July 1680 the nodal patterns on vibrating glass plates, 
running a bow along the edge of a glass plate covered with flour [47]. A hundred years later, 
the same effect was systematically studied by E. Chladni [TT]. In fact, such patterns were 
known to Leonardo da Vinci ^^, and Galileo Galilei [231 • See also some historical discussion 
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interpretation of the meaning of individual nodal deficiencies dn- It is the goal 
of this text to present one such interpretation. 

An important new approach has been developed in the last several years 
in the series of papers [I3l[311[32]. Namely, instead of concentrating on an 
eigenfunction, one can look at a partition of Vt by connected open domains 
{Dj}^^^ and try to determine whether a given partition can be the nodal 
partition of an eigenfunction of a given operator if, and if yes, what could be 
the corresponding nodal deficiency. 

Some necessary conditions on the partition are not hard to find. Indeed, one 
can introduce the graph of the partition, with each partition domain Dj 
serving as a vertex and any two nodal domains that have a (c? — l)-dimensional 
common boundary being connected by an edge. Since, due to standard unique- 
ness theorems for the Cauchy problem, the eigenfunction must change its sign 
when crossing any {d — l)-dimensional piece of the boundary of two adjacent 
nodal domains, we conclude that the following well known property holds: 

Proposition 2. The graph of the nodal partition corresponding to an eigen- 
function is bipartite. 

Another important simple observation is: 

Proposition 3. If {Dj}"-^ is the nodal partition corresponding to the eigen- 
function ipn with eigenvalue A„, then for each nodal domain Dj, one has 

(Here, as before, \k{D) denotes the kth Dirichlet eigenvalue in the domain D.) 

Indeed, by definition, ipn does not change sign in Dj and thus is proportional 
to the groundstate for H{Dj). 

This observation leads to the following notion that plays a crucial role in 
what follows: 

Definition 4. A partition P = {Dj}'^^^ is said to be an equipartition, if the 
lowest eigenvalues of all operators H{Dj) are the same, i.e. 

(2) Ai(Z}i) = Ai(Z}2) = --- = Ai(Z},). 

As we have already mentioned, every nodal partition (i.e., the partition 
corresponding to an eigenfunction) is an equipartition. 

This observation has lead to the following construction: given a natural 
number u, consider the "space" of "arbitrary" (with some natural restrictions) 
zz-partitions {Dj}, i.e. partitions with u sub-domains. Let us also introduce 
the functional 

(3) A({D,}) := max Ai(D,) 

on this space. 

One can look now at the minimal partitions, i.e. the partitions that mini- 
mize the functional A{{Dj}) for a given u. Such minimal partitions are known 
to exist [13] and the following important result holds: 



Theorem 5 (Helffer, Hoffmann-Ostenhof, Terracini [31j). Minimal bipartite 
partitions are exactly the nodal partitions of Courant sharp eig en functions. 
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Figure 1. A generic partition 

Our aim is to understand whether there is something that distinguishes 
the nodal partitions of the eigenfunctions that are not Courant sharp (which 
are the overwhelming majority) and what determines their nodal deficiencies. 
It was shown in a recent work [5] that in the quantum graph situation this 
question can be answered (see also [6] for the discussion of the discrete graph 
case). Inspired by this development, we address here the eigenfunctions of the 
operator H defined above. As we will see, one has to consider not only the 
minima, but all the critical points of the functional A on the "manifold" of 
equipartitions. Furthermore, it turns out that the Morse index of a critical 
point coincides with the nodal deficiency of the corresponding eigenfunction. 

Before we formulate our assumptions we mention that generically with re- 
spect to perturbations of the potential V{x) and/or the domain, the following 
conditions are satisfied (see [2|[3l,l36 |im - [i6| [5T] and references therein). 

Definition 6. We say the eigenfunction ipn of H is generic if 

(1) The corresponding eigenvalue A„ is simple 

(2) Zero is a regular value of the eigenfunction ipn inside fl (i.e., Vipn{x) 7^ 
whenever ipn{x) = 0). The normal derivative dtpn/dn of the eigen- 
function ipn on the boundary of Q has zero as its regular value (i.e., 
the tangential to dQ gradient of dipn/dn does not vanish whenever 
dijjn/dn{x) = 0). 

(3) The nodal set A^ = N{tpn) is the finite union of non-intersecting smooth 
hyper-surfaces 

(see Fig. [T]), where 

(4) Each Ck is a closed smooth hyper-surface in Q. 

(5) Each Bi is a smooth hyper-surface, whose boundary lies in dQ; Bi 
intersects dQ transversally. 

In fact, the intersections of Bi with dQ are orthogonal, but we will not need 
to use this information. 

We will be dealing with the generic situation only, in the sense that all 
conditions (l)-(5) above are satisfied for the eigenfunction in question. Parti- 
tions that we consider are required to satisfy similar assumptions. 
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Definition 7. A partition P = {Dj} of fi will be called generic, if its nodal 

set N = Q\ [J Dj satisfies conditions (3)- (5) above. 

We now briefly describe the constructions and results of the paper. 

Given a generic i/-partition P and a sufficiently small positive number p > 
0, we will introduce the set Vp of i/-partitions that are "close" to P in an 
appropriate sense. Here p indicates a measure of closeness. The set Vp will be 
equipped with the structure of a real Hilbert manifold by identifying it with 
a ball in an appropriate functional Hilbert space. We will denote by £p the 
subset of Vp that consists of equipartitions. 

On the set Sp of equipartitions one can consider the functional 

A :£ ^M 

that maps a partition P into the (common) lowest energy Xi{Pj) of any sub- 
domain Pj. The notation A does not contradict the one used previously in ([3]), 
since on equipartitions the two functional obviously coincide. We will also 
need some other extensions of the functional A from the set Sp of equipartitions 
to the whole Vp. Let c = {ci, . . . ,c^) G M'' be a unit simplex vector, i.e. such 
that Cj > and ^ Cj = 1. We define the functional Ac on Vp as follows: 



A,(P) = 5^c,Ai(P, 



where P EVp and Pj are the sub-domains of this partition. It is obvious that 
for any unit simplex vector c the restriction of Ac to Sp coincides with A. 
We will need the following auxiliary result: 

Proposition 8. 

(1) For any c, the functional A^ on Vp is C°° -smooth. 

(2) For a sufficiently small p, Sp is a smooth sub-manifold of Vp of co- 
dimension V — 1. 

(3) The functional A on Sp is C°°-smooth. 

This allows us now to formulate the first of the two main results of the 
paper: 

Theorem 9. Let P be a generic bipartite equipartition of a smooth domain Q. 
Then, the following statements are equivalent: 

(1) P is nodal (i.e., P is the nodal partition of an eigenfunction ip of H). 

(2) There exists a vector c = (ci, . . . ,c^) G W , Cj > 0, ^ Cj = 1, such that 
P is a critical point of the functional Ac on Vp. In this case, 

|2 



Cj = IIV^IIl2(p,) = / \'^ix)\ dx. 

(3) P is a critical point of the functional A on Sp. 

If VL is simply connected, the above statements are also equivalent to the fol- 
lowing, 

(4) At any of the boundary surfaces Cj and Bi, the normal derivatives of 
the groundstates for the two adjacent sub- domains are proportional. 
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Remark 10. If the domain Q is simply connected, the partition graph of a 
generic partition is a tree and thus bipartite automatically. Indeed, if the 
partition graph is not a tree, we can remove at least one edge while keeping 
the graph connected. That means that Q will also remain connected after 
removal of the corresponding piece of the partition boundary, Cj or Bi (call 
it B, for simplicity). However, the hyper-surface B is either closed or has a 
boundary in dQ. Therefore, due to simply-connectedness of Q, its removal 
should disconnect fi. 

Let P be the partition corresponding to a generic eigenfunction. Theorem M 
implies that, up to corrections of higher order, the functional A is given by 

A(P + AP) = A(P) + F2{AP, AP) + higher order terms, 

where F2 — the Hessian — is a quadratic form on the tangent space of Sp 
at P (which will be described in detail in Section |3]). The Morse index of 
P is then defined as the maximal dimension of the subspace on which F2 is 
negative definite. Informally, the Morse index counts directions in which the 
function A is unstable, i.e. decreasing in value. We also define the fi^ - index 
as the maximal dimension of a subspace on which F2 is non-positive. If the 
Morse and /i° indices coincide, the critical point is called non-degenerate. 
Our second main result is the following interpretation of the nodal deficiency: 

Theorem 11. Letipn be a generic (see Definition\M) eigenfunction of H and P 
he its nodal partition. Then P is a non- degenerate critical point of A restricted 
to £p and the nodal deficiency dn = n — u^^ is equal to the Morse index of A 
at the point P. 

Remark 12. To summarize Theorems [9] and [TTl instead of looking at the 
minimal points of A, one has to look at the critical min-max points, where the 
maximum is taken over a subspace of dimension equal to the nodal deficiency 
dn- This explains, in particular, why the minimal partitions correspond the 
Courant sharp eigenfunctions only. 

The structure of this article is as follows. Section [2] contains a brief expo- 
sition of the well known Rayleigh-Hadamard formula for the derivative of an 
eigenvalue with respect to the domain variation. The manifolds Vp and Sp are 
introduced in Section [31 where also the Proposition E] is proven. Theorems 
M and [TT] are proven in sections S] and El Section El contains final remarks 
and conclusions. In particular, it offers various possible generalizations of the 
results of this paper. 

2. Domain variation formulas 

In this section we provide the formulas for eigenvalue perturbation due to 
domain variation, which will be important for our considerations. Such formu- 
las have a long history, going back to J. Rayleigh [50] and J. Hadamard [29] 
and are still being developed (see [20H221 [2111251 [281 |33l|36l [391 SH] for further 
results and references). 

Let D be a proper sub-domain of Q (i.e., the closure of D belongs to Q) 
with a smooth boundary C = dD. Later on in this text D will be one of 
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Figure 2. Domain D cut out from Vthy B and its deformation. 

the sub-domains of a partition. We denote by ipi{D) and \i{D) the positive 
normahzed groundstate and the corresponding eigenvalue of H{D). 

We are interested in the variation of \i{D) with respect to infinitesimal 
smooth deformations of the boundary C . To make it precise, let us consider 
the unit external normal vector field on C and extend it into a neighbor- 
hood ?7 of C to a smooth unit length vector field N{x) whose trajectories 
are the normals to C. Let us now also have a sufficiently smooth real val- 
ued function f{x) in U . Consider the normal to C vector field f{x)N{x) and 
the corresponding evolution operators Gt of the "time" t shift along the tra- 
jectories of this field. They are defined for sufficiently small values of t and 
produce deformed surfaces Ct = GtC and the variable domains Dt bounded by 
these surfaces. Correspondingly, the ground state eigenvalue \i{Dt) is a func- 
tion of t. The following result (Rayleigh-Hadamard formula) is well known 
(e.g. [20H22l[2l[25l[28l[29l[36l[39llil[5n]) for the case when V{x) = 0. However, 
its proof (e.g. the one in [28]) is valid for non-zero potentials as well. 

Theorem 13. The t-derivative at t = of the eigenvalue Xi{Dt) is given by 
the formula 

c 

where d/dn denotes the external normal derivative on C and ipi{D) is, as 
before, the normalized Dirichlet ground-state of the domain D. 

Remark 14. If one uses a smooth vector field M{x) instead of the unit normal 
vector field N{x), an analog of formula (JH) for the t-derivative at t = of the 
eigenvalue Xi{Dt) is given by 

(5) a; = -|(^^M^^ f{x)M{x)-N{x)dS, 

c 

where d/dn denotes the external normal derivative on C and M ■ A^ is the 
inner product of vectors M and N. 

Some of the nodal sub-domains will reach the boundary. We thus also need 
to consider the case oi D (Z Q that is cut out from Q hj a. smooth surface 
(curve when d = 2) B transversal to the boundary dVt (Fig. [2]). 
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Consider again the external unit normal vector field N{x) to B. We can 
assume that this field is modified near the boundary dfl to a smooth vector 
field M{x) that is (i) non-tangential to B and (ii) tangential to dQ on the 
intersection B fl dQ. It can then be extended to a smooth vector field M near 
B such that the trajectories that start at the boundary points (i.e., points of 
B n dQ) stay on the boundary dQ and such that the level sets GfB for small 
values of t are transversal to dQ. Let / G ^^{B) (for a sufficiently large s) be 
a real valued function on B. Then, similarly to the case of an internal part of 
Vt considered above, one can deform the surface i? by a vector field /M, which 
will define sub-domains Dt of Vt with their boundaries Bt transversal to dVt. 
With this setup, the variation formula completely identical to ([5]) holds. 

Although in this case the variation formula may or may not have been 
written before, the proof of Theorem [T3] given, for instance, in P^ carries 
through without any change. 

Remark 15. If the boundary of the domain D is not connected, then the same 
variation formula holds, involving the sum of integrals over each connected 
component Cj and Bi of the boundary. 

3. Manifolds of partitions. Proof of Proposition [S] 

Let us now consider a generic //-partition P = {Pj}^^^ of Q. We need to 
introduce a manifold structure into the space of "nearby" //-partitions. The 
previous section suggests a simple way of doing so. Namely, let {Ck,Bi} be 
the smooth connected surfaces constituting the boundaries between the sub- 
domains Pj inside Q {dQ also contributes to the boundaries of some of the sub- 
domains, but is not taken into account, since it is not going to be changed). Let 
us fix a smooth non-tangential to {Ck,Bi} vector field M in a neighborhood 
U of (UfcCfc) IJ (UiBi), which satisfies the conditions imposed in the previous 
section (e.g., one can assume that outside of a neighborhood of dfi, this is 
just the unit normal field N{x), which is smoothly modified near dfl to be 
tangential to dQ). Such a field always exists under the genericity condition 
imposed on the partition P. 

Let us pick a sufficiently large positive number s (e.g., s > (rf + 4)/2 will 
suffice) and consider the space 

(6) T:= UH-'iCk)]^(®H%B, 

where -fP is the standard Sobolev space of order s. We also introduce a 
continuous linear extension operator 

(7) E : 7 ^ H'+^/^{U). 

In other words, the restriction of E{f) to C^ coincides with /. Such an exten- 
sion operator is well known to exist (see, e.g. [H]). We will also assume that 
all the extended functions E{f) vanish outside of a small neighborhood of the 
nodal set, which can be achieved by multiplication by an appropriate smooth 
cut-off function. 

Consider the ball Bp of a small radius p > around the origin in the space 
J-". Let f E Bp and Gf be the shift by time t = 1 along the trajectories 
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of the vector field E{f){x)M{x). For a sufficiently small p > 0, Gf is a 
diffeomorphism, which preserves the boundary dQ. Its action on the surfaces 
Ck and Bi leads to another zz-partition GfP of Q, which is close to the original 
partition P. We will denote this set of partitions by Vp and identify it with the 
ball Bp in the Hilbert space J-". This, in particular, introduces the structure of 
a Hilbert manifold on Vp. Notice also that we will use consistent numbering 
of the sub-domains of partitions in Vp. Namely, 

{GfP), = GjiP,), 

where, as before, for a partition 11 we denote by Uj its jth sub-domain, and 
Pj are the sub-domains of the original (unperturbed) partition P. 

We introduce now a mapping S from Vp into W^ as follows: for a partition 
n we define 

(8) s(n):=(Ai(ni),Ai(n2),...,Ai(n,)). 

Let 



A:=(A,A,...,A) 



be the diagonal in 



Definition 16. The set £p consists of all equipartitions in Vp, i.e. U & Vp 

satisfying 

Ai(ni) = Ai(n2) = --- = Ai(n,). 

To put it differently, 

(9) Sp = S-i (A) . 

One notices that the restriction of the mapping H to Sp is essentially the 
functional A of the Introduction. More precisely, 

H(n) = (A(n),...,A(n)). 

We are ready now to prove the Proposition [HI 

3.1. Proof of the Proposition [81. We prove first the following auxiliary 
result, which immediately implies the first statement of the Proposition [8] 

Lemma 17. For a small p > 0, the mapping S : Vp ^M." is C°°. 

Proof. It is sufficient to prove that for any j, the mapping 

/ ^ XiiGfiP,)) 

is smooth as a mapping from the ball Bp to M, if p is sufficiently small. Thus, 
one can restrict attention to a single sub-domain D. The often employed in 
such circumstances idea is to replace domain dependence with varying the 
coefficients of the differential operator in a fixed domain. Then the smooth 
dependence of Ai becomes a standard perturbation theory result (e.g., [25|[3l]). 
Let us consider first the case of a sub-domain D that does not touch the 
boundary. Consider the mapping 

^f : X -)■ y:= x + E{f){x)M{x) 



10 GREGORY BERKOLAIKO\ PETER KUCHMENT\ AND UZY SMILANSKY^-^ 

of the domain D into Vt. For sufficiently small p > 0, it is a C^-diffeomorphism 
of D onto a "nearby" sub-domain D* of Vt. The quadratic form of the operator 
H{D*) is given as 



du 



dy 



+ V{y)\u{y)\'\dy. 



D* 

Changing variables back from y to x, we arrive at an operator Hf{D) in the 
fixed spatial domain D, but with variable coefficients now: 

{Hfu){x) = V ■ Af{x)Vu{x) + Vf{x)u{x), 

where the matrix- valued function Af{x) is of class C^, / — )■ A/ is a C°°- 
mapping from H'^ to the space of C^-matrix functions, and / — )■ V/ is a C°°- 
mapping from H'^ to C^. 

We have thus replaced the domain dependence with the smooth dependence 
of the coefficients of the operator. The operator Hf acts continuously from 
H'^{Q) to L^(fi) and for / = coincides with H{D). Moreover, the mapping 
Af — )■ Hf is a continuous linear mapping from the space of C^-matrix func- 
tions to the space of bounded operators from H'^{Q) to L^(f2). Thus, for a 
sufficiently small p we get a smooth family of Fredholm operators between the 
aforementioned spaces. Due to the simplicity of \i{D), the standard perturba- 
tion theory shows that Xi{Hf) depends smoothly on /, for a sufficiently small 
radius p. 

A similar consideration works when D reaches the boundary of Q, i.e. at 
least one of the boundaries Bi is involved. Without loss of generality, we can 
assume that only one such Bi is involved. Introducing an appropriate smooth 
coordinate change, one can reduce consideration to the cylinder B x [—e,6) 
for small e > 0, with dB x [~e,e) as the corresponding part of dQ. Then the 
same reduction to a fixed domain but varying operator as before is possible, 
which again implies smooth dependence of Ai on /. D 

Let us now address the second statement of the Proposition, that Sp is a 
smooth sub-manifold of Vp of co-dimension z/ — 1. We will employ for this 
purpose the formula iQ in conjunction with the domain variation formula, 
equation (|5]), and a transversality theorem. 

According to Lemma [T71 the mapping S : Pp — )■ M*^ is a smooth mapping of 
Banach manifolds. The pre-image of the diagonal A C M'^ coincides, according 
to 0, with Sp. We would like to know whether this pre-image is a smooth 
sub-manifold, and of what co-dimension. This is exactly the question tackled 
by the transversality theorems. Namely, if we can show that the mapping S 
is transversal to the diagonal one-dimensional sub- manifold A of M'^, this will 
prove that the pre-image of A is a smooth sub-manifold of co-dimension z/ — 1. 

Definition 18. (e.g., [T||9|B0]) The mapping E : Vp ^ W is transversal to 
A, if at any point ( = S(t>), v G Vp, that belongs to A, the vector sum of the 
tangent space T(^A to A at C and of the range D'E{TyVp) of the differential D'E 
on the tangent space T^Vp is the whole space W . 

Theorem 19. (e.g., |21 Sect. 3, Theorem 2] or f^ Sect. 5.11.7]) 

IfS is transversal to A, then Sp = S~^(A) is a smooth sub-manifold ofVp of 

co-dimension z/ — 1. 
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Thus, to finish the proof of the second statement of Proposition [8], it only 
remains to prove the transversahty of the mapping S to the diagonal A. The 
Rayleigh-Hadamard domain variation formulas are helpful here. 

The tangent space to the diagonal is spanned by the vector (1, 1, ... , 1). We 
will demonstrate that the range of -DS contains a subspace of dimension at 
least V — 1. If the dimension is z/, there is nothing further to prove, but if the 
dimension is i/ — 1 we will show that the vector (1, 1, ... , 1) does not belong 
to the subspace. 

Consider the partition graph F that corresponds to an (automatically generic) 
zz-partition P E Vp. The vertices of the graph correspond to the sub-domains 
Pj and the edges to the interfaces Cj,Bi. We will identify the target space 
W of the mapping S with the space of real valued functions on the set V of 
vertices of the graph F. Consider a pair of adjacent sub-domains Pi and Pj 
with the common part of their boundary 5* (one of Cj,Bi). We consider the 
corresponding vertices Vi,Vj and the edge s of F. We restrict our attention 
to functions f E J-' that are non-zero on S only and find the corresponding 
directional (Gateaux) derivative of H at P in the direction of /. Equation ([5]) 
shows that the only non-zero components of this derivative correspond to ver- 
tices Vi and Vj. Since M(x) is non-tangential to S, M{x) ■ N{x) is sign-definite. 
Choosing / of the same sign we get 

Df'Ei = - j (^^p-\ f{x)M{x)-N{x)dx < 

and, similarly, DfEj > 0. Here we assumed that the normal N{x) is directed 
outward with respect to Pi and, correspondingly, inward with respect to Pj. We 
also used the fact that, due to the standard uniqueness theorems, dipi{Pi) / dn 
is not everywhere zero. 

Repeating this procedure for every pair of adjacent sub-domains we arrive 
at a collection of vectors, one for each edge of the partition graph, that have 
two non-zero components of the opposite sign each. To characterize the space 
spanned by these vectors, we arrange them as rows of a matrix and find its 
kernel. Due to the connectedness of the graph, the kernel is at most one- 
dimensional. If the kernel is empty, the vectors we found span all of W^. If 
the kernel is spanned by a vector u (this is the case if the domain Q is simply 
connected, as this implies that the partition graph is a tree), then u must have 
all components of the same sign. The vector u is therefore not orthogonal 
to the vector (1, !,...,!), and the latter vector complements the derivative 
vectors to span W. 

This finishes the proof of transversahty and thus of the second statement of 
the proposition. Since the first two claims of the proposition imply the third 
one, the proof is completed. 

4. Proof of Theorem M 

4.1. Proof of the equivalence ([I]) <^ ({2]). If P is the nodal partition of a 
real- valued eigenfunction ip, then, as we have already mentioned before (Propo- 
sition |3]), the restrictions of ip to the nodal domains are proportional to the 
groundstates in these domains. Denote these proportionality constants by 
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aj. Since the eigenfunction ip is continuously differentiable, the groundstates 
ipilPj) scaled with the corresponding factors aj have matching normal deriva- 
tives at the common boundaries: 



(10) a, 



dn 



(Jj'j 



dn 



where sub-domains Pj and Pj have the common boundary S, and n is a normal 
vector to S. 

Now let Cfc = a1 and consider the Gateaux derivative of the functional Ac 
in the direction / that is non-zero only on a single boundary between the 
sub-domains Pi and Pj. Since the only affected terms in Ac are afX{Pi) and 
OjX^Pj), the derivative is 



(11) aU ^-^^^^ -aU --^^^^ f{x)M{x) ■ N{x)dS^ 




where we applied the Rayleigh-Hadamard formula, equation ([5]). The differ- 
ence in signs arises since N{x) points outward with respect to Pi but inward 
with respect to Pj. Now we observe that, due to ( TTOl) . the integrand is iden- 
tically equal to zero and thus the Gateaux derivative in the direction of / is 
equal to zero. The same is obviously true for arbitrary variations /, involving 
any number of boundaries. Thus, the nodal partition is a critical point of the 
functional Ac. 

Conversely, if a partition Pp is a critical point of Ac, we get that the Gateaux 
derivative of Ac is zero in any direction f{x). This implies the equality 



dn J \ dn 

on the common boundary S of any two neighboring domains Pi and Pj . Setting 
ttfc = '^\fck and choosing the signs so that any two neighboring domains have 
different signs (possible due to bipartiteness) ensures that (TTUl) is satisfied. 
Then the function ■?/' defined by 

V" |Pfc= akipi{Pk) 
is an eigenstate of H. 

4.2. Proof of the equivalence ([2]) <^ (|3]). If P is a critical point on Vp of 
the functional Ac, then the restriction of Ac to Sp is a critical point on Sp. But 
on Sp any functional Ac coincides with A. 

Conversely, assume that P is a critical point of A on Sp. We can extend 
the functional A to the whole Vp as Ai(Pi). Since Sp can be given by the 
smooth relations Ai(Pi) — Ai(Pj) = for j = 2, 3, . . . , z/, the Lagrange mul- 
tiplier method implies that P must be a critical point of a non-trivial linear 
combination A^ := Yl^j-^ii^j)- ^11 ^j ^^^ of ^^e same sign: otherwise there 
are two neighboring domains with bj of different signs and the variation of A^ 
with respect to the boundary between the two domains cannot be zero for a 
sign-definite /(x) (see equation (ITT]) ). Thus the vector of coefficients bj can be 
normalized to be a unit simplex vector. This finishes the proof of Theorem [91 
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5. Proof of Theorem [TT] 

Let us present first the strategy of tlie proof. As we have already seen, it 
is sometimes useful to play with different extensions of the functional A from 
the (local patch of the) space £p of equipartitions to a larger manifold. While 
previously it was the (local patch of the) space of all partitions, now we need 
some further extension. Indeed, we would like to compare somehow the nodal 
count V (which is fixed) with the consecutive number n of an eigenfunction ipn- 
It is hard to observe where the information about n is hidden in the spaces 
of partitions themselves. On the other hand, the quadratic form (or Rayleigh 
quotient) contains this information. To use it we will extend the functional A 
to a larger space that is a functional space, not just a set of domain partitions. 
Then we will have to restrict back in order to compare the Morse indices of A 
and of its extension. 

Before implementing this program, in the following sub-sections we start 
proving some auxiliary statements that will come handy later on. 

5.1. Critical points on direct sums of spaces. 

Theorem 20. Let X = Y ^Y' be a direct decomposition of a Banach space. 
Let also f : X ^ W be a smooth functional such that (0, 0) G X is its critical 
point of Morse index m. 

If for any y in a neighborhood of zero in Y , the point {y, 0) is a critical point 
of f over the affine subspace {y} x Y' , then the Hessian F2 of f at the origin, 
as a quadratic form in X, is reduced by the decomposition X = Y ^Y' . 

In particular, let Y be the locus of minima of f over the affine subspaces 
{y} X Y', I.e. 

(12) (y, 0) = arg min f{y, y'), 

y'eY' 

for any y in a neighborhood of zero in Y . Then the Morse index of as a 
critical point of the restriction /|y is equal to m (i.e., the same as the Morse 
index of this point on the whole neighborhood of zero in X). 

Finally, if (0, 0) is a non- degenerate critical point of f on X , then is 
non-degenerate as a critical point of fly. 

Proof. We can assume, without loss of generality, that /(O, 0) = 0. Using this 
and the condition of the criticality of the origin, the Taylor formula of the 
second order for / on X near the origin is 

/(y, y') = My, y) + B{y', y') + C{y, y') + higher order terms. 

Here A is a quadratic form in Y , B is a, quadratic form in Y' and C is a bilinear 
form acting on F x Y' . We will now take the Gateaux derivative of /(y, y') in 
the direction z' G Y' and evaluate it at (y, 0). 

The derivative of A{y,y) is zero (since A does not depend on y'). The 
derivative of B{y',y') is B{z',y') + B{y',z'), which vanishes after substituting 
y' = 0. Thus we find 

D^o,.')fiy,0) = Ciy,z'), 

which must be zero since {y, 0) is a critical point over {y} x Y'. Since z' and 
y are arbitrary (provided y is sufficiently close to 0), the bilinear form C is 
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identically zero. The Hessian is thus 

F^ = A{y,y) + B{y',y'). 

This proves the first statement of the Theorem. 

If condition (TT2l) is fulfilled, then B{y',y') is non-negative definite and the 
negative subspaces of F2 and A coincide (after the natural projection). 

Finally, if (0,0) is non-degenerate, it implies that /i°(A) = /u(^), where /i 
denotes the Morse index. Indeed, otherwise /U°(F2) > /i°(A) > /i(^) = /w(-F2) 
which contradicts the non- degeneracy of (0, 0). D 

Remark 21. It is obvious from the proof of the theorem that its second 
statement can be generalized in the following manner: if Y is the locus of 
critical points of / of index m' (locally independent of y) then the Morse 
index of the restriction /|y is equal to m — m'. We will only need m' = in 
the present manuscript but the more general version becomes useful in other 
contexts, such as [7]. 

5.2. Some objects needed for the proof. Due to the local character of 
Theorem [m all the constructions below are needed only locally, near a generic 
eigenfunction ipn indicated in the statement of the theorem, and correspond- 
ingly near its nodal partition. The basic notions and facts concerning finite- or 
infinite- dimensional vector-bundles that we use below can be found in many 
standard sources on topology (e.g., mSS]) or in the survey |52j, where such 
bundles are studied in relation to the operator theory. 

We will be considering again the (local) manifold Vp of partitions "close" to 
the nodal partition P{ipn) and its sub-manifold Sp of codimension u — 1 that 
consists of equipartitions only. 

Definition 22. In the trivial bundle Vp x Hq{Q) over Vp, we consider a fibered 
subset B that has the fiber over a partition P consisting only of functions 
vanishing on the partition's interfaces A^ U dQ. In other words, this fiber is 

Lemma 23. B is a smooth locally trivial vector sub-bundle of the trivial bundle 

VpxH^{n)^Vp. 

Proof. The proof follows the same line as the one of Lemma [T71 Namely, the 
dependence on the partition P is replaced, using a smooth family of diffeomor- 
phisms, with a fixed partition 11. The corresponding change of variables in the 
functions from Hp maps them to H^. Then B becomes just the trivial sub- 
bundle Vp X ( . Hq{IIj) ] inVpX -f^g (^) ^ ^p- Inverting the diffeomorphisms 
provides a smooth trivialization of B, which proves the lemma. D 

Definition 24. We denote by SB the locally-trivial bundle of the unit (in 
L2-norm) spheres of the fibers of B. 

The restrictions of B and SB to Sp (clearly locally-trivial bundles) will be 
denoted by Be and SBe correspondingly. 

We will now restrict the bundle B further. 
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Definition 25. We denote by C the vector bundle whose fiber over a partition 
P consists of functions of the form ^,- Cj'?/'i(Pj), where Cj are real constants 
and V'i(-Pi) is the normalized positive groundstate on the sub-domain Pj. 

Correspondingly, SC consists of the unit (in L2-norm) spheres of C and Ce 
and SCe are restrictions of the corresponding fibered sets to £p. 

The following lemma shows that C and SC are locally-trivial bundles, and, 
in particular, C is a locally trivial vector bundle. 

Lemma 26. C is a smooth v -dimensional vector sub-bundle of B, and thus of 
the trivial bundle 

VpXHl{Sl)^Vp. 

Proof. The proof follows the same line as in Lemmas [T7] and [231 Namely, 
after applying a smooth family of diffeomorphisms, one deals with a fixed 
partition 11, but instead with the operator whose coefficients depend smoothly 
on P. Perturbation theory shows that the corresponding ground-state fj in 
each sub-domain IIj depends smoothly on P, as a vector in Hlijlj). We extend 
it, without changing the notation, by zero to the whole domain Vt. Then 
{fj} is a smoothly dependent on P frame of v linearly independent vectors 
in Hq{VL). Thus, this frame spans a smooth finite-dimensional vector-bundle. 
After applying the inverses of the diffeomorphisms, we get the claim of the 
lemma. D 

5.3. The functional A and the quadratic form of H . Consider the qua- 
dratic form on HI{VL) 

Q[f]--= [ {\Vf{x)\' + V{x)\fix)\')dx. 

It can, by restricting to each fiber, be defined as a smooth functional on the 
vector bundle B and its sub-bundles that we considered above. 

Lemma 27. The point {P{ipn),4'n) ^^ SB is a non- degenerate critical point 
of Q of Morse index /i = n — 1. 

Proof. It is clear that Q on Hq{Q) has an n — 1-dimensional subspace on which 
its Hessian at ipn is negative. Namely, this is the subspace generated by the 
eigenfunctions ipi, ...,ipn-i- If we show that these directions are among the 
tangential ones to SB at (P('?/'„), ■?/'„), this will prove that /i > n — 1. 

Due to the locally-trivial structure of SB, there are two main ways to get 
tangential vectors to SB. One is to vary the partition P (which will give 
"horizontal" tangent vectors to the bundle SB). The other is to keep the 
partition fixed, while varying the function ipn, keeping the nodal set fixed 
("vertical" tangent vectors). The vertical tangent vectors are just arbitrary 
functions in Hq{Q) that vanish on the nodal set Z of the partition P{'4)n) (we 
have previously denoted this space Hp). The horizontal tangential vectors 
look as follows: 

E{f)M'4,n = E{f)Ni,n{x){M{x) ■ N{x)), 

where / G H^ on the nodal set of ipn is the function defining this set's infin- 
itesimal variation. Also, the notation Xg for a vector field X and a function 
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g (e.g., Mipn and Nipn) mean the derivative of the function g along the field 
X. Notice that, since M{x) is non-tangential to Z, M{x) ■ N{x) is a smooth 
separated from zero function on the nodal set. 

We now show that any ipj with j < n can be represented as a sum of vertical 
and horizontal vectors. We notice that the genericity condition requires in 
particular that zero is a regular value of the normal derivative of ipn on dQ. 
This implies that the derivative Nipn is a smooth function that has a non- 
degenerate zero at dQ on the nodal set Z. Since ipj is a smooth function 
vanishing on dfl, the function 



/(^) 



ijj{x) 



{NtPn)M{x) ■ N{x) 

belongs to i7* on the nodal set. Hence, the horizontal tangent vector in the 
direction of / 

h := {N^n)E{f)M{x) ■ N{x) 

coincides with ipj on the nodal set. This means that the difference g := ipj — h 
belongs to Hp, and thus is a vertical tangent vector. This shows that each 
eigenfunction ipi, ...,iljn_i can be represented as the sum of a vertical and hor- 
izontal vectors and thus is tangent to SB. This proves the estimate fi > n — 1 
for SB. 

We will now prove that the index /i° cannot exceed n — 1. This will verify 
that yU = n — 1 and that the critical point is non-degenerate. 

Suppose that there is an n-dimensional subspace L in the tangent space 
to SB at {P{%l)n))'4'n)) where the Hessian of Q is non-positive. Since each 
fiber of SB consists of functions from the space Hq{VL)., there is a tautological 
mapping {P, f) i— ?■ / from SB into Hq{Q) (in fact, into the set of functions 
of unit L2-norm). If the Frechet derivative of the tautological mapping has 
zero kernel, then the subspace L will produce an n-dimensional subspace of 
Hq[Q) transversal to ipn, where the Hessian of Q at ipn is non-positive, which 
is a contradiction. So, let us show that the kernel of the Frechet derivative is 
zero. Due to the local trivial structure of SB, one sees that the image of any 
tangent vector under the Frechet derivative has the form 

g = E{f){x){M^Pn)+v{x), 

where the function / G J-" is responsible for the infinitesimal variation of the 
nodal set of ipn, E is the previously introduced extension operator from the 
nodal set to Q, and a function v G Hp{Q) corresponds to the infinitesimal 
variation in the fiber direction (i.e., the pair (/, f) describes a tangent vector 
to SB). Suppose now that g = 0. In particular, g\z = 0. Taking into the 
account that v\z = (which is true for any function from the space Hp), we 
conclude that 

E{f){x) (M^„) \z = f{x)M{x) ■ VMx)\z = 0. 

Since Vipn{x)\z is nonzero and orthogonal to Z, and M is transversal to Z, 
we see that 

M{x)-VMx)\z7^0. 

This implies that f{x) (defined on Z only) vanishes identically. Due to the 
linearity of the extension operator, the extension E[f) vanishes everywhere in 
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Q. Hence, g = v. Since g = 0, we conclude that v = and thus the Frechet 
derivative of the tautological mapping is injective. This finishes the proof of 
the lemma. D 

Lemma 28. The point {P (ipn) , 'ipn) in SC is a non- degenerate critical point 
of Q of Morse index n — 1. 

Proof. We know that the eigenfunction ipni^) has the form 

V 

i=i 

where {Dj\ is the nodal partition of ■^/'n and each ground state ipi{Dj) is 
extended by zero to the whole domain Q, X] c? = 1, and none of the coefficients 
Cj are equal to zero. Thus the coefficients Cj remain bounded away from zero 
in a neighborhood U of i/jn in SC. 

As we have already seen, C is a smooth z/-dimensional sub-bundle of B and 
thus SC is a smooth sub-manifold of SB. We will now introduce near the 
point {P{ipn),'4^n) a smooth foliation of SB by manifolds transversal to SC, 
such that SC will be the locus of minima of Q over the leaves of this foliation. 
Then, after a smooth local change of coordinates, we will be in the situation of 
Theorem l20l Thus, the Morse index of Q at {P{ipn),i^n) on SC will be equal 
to the one on SB, which will prove the first claim of the lemma. 

So, let us finish the proof by constructing such a foliation. Consider the 
following mapping T from a neighborhood U of (P(^„), ipn) ^ SB to Vp x S'''~^, 
where S'^~^ is the unit sphere in M.'^: 

T(P,/) = (P,||/|U,(P,),...,||/|U,(P^)). 

Notice that none of the components ||/||l2(-P) of ^^e vector T(P, /) vanishes 
(since this is the case for T(P(-?/)„), ■?/'„)). 

By the arguments provided before, this is a smooth mapping. It is also clear 
that it is a submersion. 

Let w = (P, X]cj'?/'i(Pj)) be a point in U H SC near {P{ipn),ipn)- Consider 
the leaf L^ = T~^(P, ci,C2, . . . ,0^). Due to the submersion property of T, 
the leaves Ly^ form near {P{ipn),ipn) a smooth fibration. Since the differen- 
tial of T on SC is surjective, this foliation is transversal to SC. Moreover, 
the groundstate ipi{Pj) is the location of the minimum of the corresponding 
Rayleigh quotient, and therefore 

Q\f\ = Y.Q{f\p;^>Y.Q\c,^[p,)\, c, = ii/iu,(p^,). 

i=i i=i 

That is, the minimal value of Q on the leaf L^ is attained exactly when /,■ = 
Cjipi^Pj) for all j, i.e. on SC. D 

Lemma 29. The point {P (ipn) , "ipn) is a (degenerate) critical point of Q in 
SCe with the Morse index not less than n — v and the ^^ -index not more than 
n — 1. 

The restriction of the quadratic form Q to the fibers of SCe is the pull-hack 
of the functional A from the base Ep. In other words, if vr : SCe ^~^ £p is the 
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bundle projection, then for any x G SCe, 

Q{x) = ^{'K{x)). 

Proof. According to Proposition |H1 SCe is a smooth sub-manifold of SC of 
co-dimension v — 1. Hence, upon restriction to SCe-, an index (either Morse 
or /i°) cannot decrease by more than z/ — 1. A restriction also can not increase 
an index. These observations, together with Lemma [281 prove the first two 
claims. 

Let / = Xlj '^jV'i(-Pj) be an element of the fiber of SCe over a partition 
P = {p-}. Then 

Since P is an equipartition, all the values Xi{Pj) are equal to the same value 
A(P). Taking into account that X] kjP = ll/llL2(n) ~ -'-' ^^ S^^ 

which proves the last statement of the lemma. D 

5.4. Index of A and nodal deficiency. We are now ready to prove Theorem 
[TTlby obtaining estimates for the Morse index /x from two sides. 

Proof of the estimate from below: fi > n — z/^„. From Lemma [29] we see that 
the pull-back of A to SCe has Morse index at least n — u at the point 
{P{ilJn),4'n), i-e. there is a tangent subspace of dimension at least (n — z/), 
where the Hessian of Q is strictly negative definite. On the other hand, the 
pullback is constant along the fibers of vr : SCe ^ ^p- Therefore the above 
subspace must correspond to such a subspace for the Hessian of A. D 

Proof of the estimate from above: ^ < fJ^^ < n — u^^^. We use Lemma [21] again. 
It shows that the index fi^ of A at P{ipn) cannot exceed (n — 1) — (z/ — 1) = n — u. 
Indeed, the pullback of A to SCe add a (z/ — l)-dimensional subspace on which 
the Hessian is zero. But the fi^-index on SCe is bounded hj n — 1. D 

To summarize, together with the trivial inequality /x < /i° (see the discussion 
prior to Theorem [TT]) we have shown that 

This shows that the critical point is non-degenerate (i.e. /i = /i°) and of Morse 
index fi = n — z/^„ , finishing the proof of Theorem [11] 

6. Final remarks and conclusions 

(1) The results of this paper (Theorems [5] and [TT]) translate without any 
changes in their proofs to the case when fi is a compact smooth Rie- 
mannian manifold with or without boundary. 

(2) Smoothness conditions imposed on the domain, potential, and partition 
interfaces, can certainly be weakened. In this text, we have decided not 
to do so, in order not to complicate considerations unnecessarily. 
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(3) The sets V and £ that we considered involved only generic partitions, 
which allows perturbing the boundaries of sub-domains independently, 
simplifying the structure of the manifold of partitions and the con- 
sequent considerations of the text. In general, however, the smooth 
pieces of partition manifolds are joined into singular "varieties" V and 
£, where the junctions occur when the partition interfaces start meet- 
ing each other. It would be interesting to see whether one could prove 
an analog of Theorems [9] and [11] for such non-generic partitions. The 
authors believe that something of this nature can be done. 

(4) An interesting connection between the zeros of eigenf unctions and sta- 
bility index with respect to another suitably defined perturbation has 
emerged recently. In the preprint [7] (see also [12] for an alternative 
proof), it was shown for the graph case that the "surplus" number of 
zeros (pn — {n — 1) is equal to the Morse index of the eigenvalue consid- 
ered as a function of magnetic perturbation of the Hamiltonian. Here 
4>n is the number of zeros of the n-th eigenfunction. 

(5) The reader probably has noticed that there is a significant flexibility 
in our choice of the functional (e.g., Sobolev) spaces. Probably, one 
could achieve the same results working in the C^ spaces of smooth 
functions. Our choice of the Sobolev scale was due to a known simpler 
interpretation of the Morse theory on Hilbert, rather than Banach, 
manifolds. 
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